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I.  INTRODUCTION 

In  the  early  1970s,  stealth  technology  was  revealed  to  the  public.  Since  then,  the 
term  stealth  has  generally  been  interpreted  to  mean  "invisible  to  radar."  In  reality,  however, 
stealthy  targets  are  not  completely  invisible.  Stealth  technology  merely  reduces  the  radar 
cross  section  (RCS)  of  the  target.  The  RCS  of  a  target  is  denoted  by  ct  and  may  be  defined 
as: 

Power  reflected  to  receiver  per  unit  solid  angle 

C7  = — 

Incident  power  density  /  4n 
A  great  deal  of  effort  is  spent  to  reduce  RCS  using  four  basic  approaches: 

•  Target  Shaping 

•  Surface  Material  Selection  and  Coatings 

•  Active  Cancellation 

•  Passive  Cancellation 

Each  of  these  methods  involves  trade-offs.  For  instance,  the  optimum  target  shape  to 
achieve  reduced  RCS  may  result  in  aerodynamic  problems.  Additionally,  RCS  reduction 
methods  (except  for  shaping  )  are  usually  very  narrowband  and  are  effective  over  limited 
spatial  regions.  These  factors  dictate  that  RCS  reduction  methods  must  be  chosen  based  on 
the  platform's  mission  and  potential  threats.  The  main  goal  of  RCS  reduction  is  to  achieve 
an  RCS  value  below  an  acceptable  threshold  over  a  specified  range  of  frequencies  and 
angles. 

The  relationship  between  the  electrical  characteristics  of  the  material  employed  and 
the  target's  resulting  scattered  field  is  a  very  important  concept  in  RCS  reduction  (RCSR). 
The  goal  of  materials  selection  is  to  choose  a  material  that  yields  the  smallest  reflection  of 
radar  waves.  The  use  of  radar  absorbing  materials  (RAM)  is  crucial  to  RCSR,  but  because 


of  the  complexity  of  the  problem,  RCS  reduction  using  RAM  has  historically  been  based  on 
intuition  and  experience. 

This  thesis  examines  a  method  of  synthesizing  a  target's  surface  material's  electrical 
properties  to  obtain  a  specified  RCS.  Synthesis  is  a  special  case  of  the  more  general 
problem  of  inverse  scattering,  whereby  information  on  an  unknown  target  is  derived  from 
its  scattered  field.  In  RCS  synthesis,  the  target  shape  is  usually  known;  determining  the 
surface  material  is  the  primary  goal.  The  synthesis  problem  can  be  solved  in  either  the  time 
or  frequency  domains.  This  thesis  addresses  the  problem  in  the  frequency  domain. 

The  solution  of  integral  equations  provides  a  basis  for  obtaining  surface  impedance 
information  from  a  target's  RCS.  In  general,  both  electric  and  magnetic  currents  can  exist 
on  the  surface.  Integral  equations  can  be  solved  by  using  the  method  of  moments  (MM), 
which  is  a  frequency  domain  technique.  The  MM  solution  employed  here  requires  that  the 
entire  structure  to  be  modeled  be  decomposed  into  a  number  of  surface  patches 
(subdomains),  which  must  be  small  compared  to  the  wavelength.  There  are  two  approaches 
to  deriving  the  integral  equations  for  surface  impedance  synthesis.  These  are  the 
approximate  method  and  the  rigorous  method.  In  the  approximate  method,  the  unknown 
surface  impedance  is  expanded  in  a  series  of  basis  functions  with  unknown  coefficients,  and 
an  estimate  of  the  surface  current  is  provided.  The  rigorous  method  requires  that  the 
coefficients  of  series  expansions  for  both  the  impedance  and  current  be  determined 
simultaneously. 

Previous  works  on  RCS  synthesis  have  discussed  the  problem  in  terms  of  the 
approximate  method  [1  and  2].  The  rigorous  method  is  the  focus  of  this  thesis.  In  general, 
for  the  rigorous  method,  a  nonlinear  set  of  simultaneous  equations  must  be  solved.  The 
present  work  will  be  limited  to  a  resistive  surface  and  the  target  geometry  to  a  flat 
rectangular  plate.  Additional  boundary  conditions  will  be  utilized  to  solve  the  synthesis 
equations  for  both  sets  of  unknown  coefficients.  In  this  limited  case  the  resulting  system  of 
equations  is  linear. 


Chapter  II  presents  theoretical  background  required  in  the  formulation  and 
evaluation  of  the  problem.  Reflection  and  transmission  coefficients,  boundary  conditions, 
and  surface  impedance  concepts  are  explained.  In  Chapter  III,  the  RCS  of  a  rectangular 
resistive  sheet  is  calculated  with  MM  and  the  RCS  synthesis  procedure  for  the  rigorous 
method  is  developed.  Appropriate  equations  are  derived  and  expressed  in  matrix  form  for 
the  bistatic  case.  In  Chapter  IV,  computer  implementations  of  these  synthesis  equations  are 
discussed.  Chapter  V  summarizes  the  results  of  the  synthesis  and  presents  conclusions  and 
recommendations  for  future  research. 


n.  THEORETICAL  BACKGROUND 

A.         SURFACE  RESISTIVITY 

The  goal  of  radar  cross  section  reduction  is  to  eliminate,  or  at  least  attenuate 
significant  scattering  sources  on  the  target.  Thin  layers  of  lossy  materials  are  of  special 
interest  in  RCSR  design,  primarily  in  the  treatment  of  surface  waves.  Discontinuity 
boundary  conditions  are  used  to  mathematically  represent  thin  films.  They  are  the  resistive 
sheet  and  conductive  sheet  boundary  conditions.  The  boundary  conditions  define  the 
relationship  between  the  current  and  the  tangential  field  components  on  the  surface.  A 
resistive  sheet,  which  is  the  primary  focus  of  this  thesis,  is  an  infinitely  thin  imperfect 
electric  conductor.  Its  conductivity  is  finite,  which  means  a  e  <  oo  ,  and  it  does  not  support  a 

magnetic  current  ( Jms  =  0  ).  For  resistivity,  the  unit  ohms  per  square  ( Q  /  sq )  is  generally 
used.  This  unit  is  derived  from  the  basic  definition  of  resistance  for  the  special  case  of 
length  equal  to  width  ( I  =  w ) 

£  1 

GctW        (Jct 

where  ac  is  the  conductivity  of  the  material  per  meter  and  t  is  the  thickness  of  the 

material. 

There  are  two  limiting  values  for  surface  resistance.  The  first  is  Rs  =0,  which  represents  a 
perfect  electric  conductor.  The  second  is  R=co ,  which  represents  a  surface  matched  to  free 
space. 

The  boundary  conditions  may  be  expressed  in  general  vector  form  as  [1] 

nxE(+)-nxE(-)  =  0  ,  (2) 

hxH(+)-hxH{-)  =  Js   ,  (3) 

with: 


nx{hxE(±)}=-RsJs,  (4) 

where  h  is  the  unit  vector  normal  to  the  sheet,  and  is  directed  to  the  upper  side  of  the  sheet, 

and  J5is  the  total  electric  current  supported.  Plus  and  minus  signs  represent  upper  and 
lower  faces  of  the  sheet,  respectively. 

The  electromagnetic  dual  of  the  resistive  sheet  is  the  magnetically  conductive  sheet, 
which  supports  only  magnetic  currents.  A  combined  sheet  consisting  of  resistive  and 
conductive  sheets  is  equivalent  to  an  impedance  surface,  which  is  a  surface  where 
permittivity  and  permeability  both  differ  from  the  surrounding  medium  [1].  The  surface 
impedance  approximation,  also  known  as  the  Leontovich  boundary  condition,  is: 

^iK  (5) 

where  E^  and  H^  are  the  total  tangential  fields  on  the  surface,  and  tjs  is  the  surface 
impedance  in  ohms  [3].  Solution  of  a  scattering  problem  for  the  resistive  sheet  gives  the 
solution  for  the  conductive  case  by  way  of  the  duality  principle. 

B.         REFLECTION  COEFFICIENT 

When  a  traveling  wave  reaches  an  interface  between  two  different  regions,  it  will  be 
reflected  or  transmitted  with  the  relative  magnitudes  of  the  two  components  determined  by 
the  impedances  of  the  two  regions.  The  reflection  coefficient  for  a  resistive  thin  film  can  be 
derived  from  the  boundary  conditions  in  equations  (2),  (3),  and  (4).  The  geometry  is  shown 
in  Figure  1 .  Bistatic  scattering  for  parallel  and  perpendicular  polarizations  are  investigated 
separately.  Specular  reflection  can  be  assumed  because  the  sheet  is  infinitely  thin. 

1.  Parallel  Polarization 

For  parallel  polarization  (TM  case),  the  incident  wave  is 

E,  =  &#«-*' ,  (6) 

where  k  is  the  incident  wave  vector.  After  some  mathematical  manipulations,  the  reflection 


coefficient  is  found  to  be:  [1] 


r     = 

1  TM 


-x\0  cos0;. 


2^5+r|0cose,  ' 
The  surrounding  medium  is  assumed  to  be  free  space  with  impedance  r|  0 . 


(7) 


Figure  1 .  Reflection  and  Transmission  from  a  Resistive  Sheet 
2.  Perpendicular  Polarization 

For  perpendicular  polarization  (TE  case)  the  incident  wave  is: 

and  the  reflection  coefficient  for  the  resistive  sheet  is    [1] 


(8) 


r     = 

1  TE 


2^cos6,. +Tj0 


(9) 


C.         RADIATION  INTEGRALS 

The  radiation  integrals  are  integral  solutions  to  Maxwell's  equations,  and  are  also 
called  Stratton-Chu  integrals.  They  can  be  derived  directly  by  taking  the  curl  of 
Maxwell's  first  two  equations,  using  the  vector  of  Green's  theorem,  and  then  integrating 
[3].  E-  and  H-fields  are  determined  from  current  distributions  by  using  the  radiation 
integrals.  On  a  resistive  sheet  the  magnetic  currents  are  zero.  The  field  equations  are 
further  simplified  when  only  the  far  zone  fields  are  of  interest. 

In  the  far  zone  ( kr  »  1 )  and  the  E-  and  H-field  vectors  are  orthogonal  to  each 
other  and  the  direction  of  propagation  resulting  in  TEM  mode  fields.  The  equations  used 
in  this  section  are  derived  by  Balanis  [4].  A  rectangular  coordinate  system  is  used  in 
formulating  the  solution  for  the  rectangular  sheet  problem.  The  geometry  is  shown  in 
Figure  2. 


Figure  2.  Coordinate  System  for  Scattering  from  a  Planar  Rectangular  Plate 


The  far-field  spherical  components  of  the  E-field  are: 


and 


where 


jke~jkr 
£e=— ^V^e,  (10) 


J7  Jke~Jkr       AT 

^~~~^r^°  *  (11) 


NQ  =  U^cosecos^  +  ^cosesirK))]  e^'^dx'dy' ,        (12) 
N*  =  J{  [~  J*  sin+  +  Jy  <*»♦]  ejkr,cosxvdx'dy'    ,  (13) 


and 


r'cosij/  =  x'sin0cos(})  +  7'sin9sin(j)  =  x'w  + >>'v  =  g.         (14) 
The  direction  cosines  are 

w  =  sin0cos(|)  , 

v  =  sinGsincj)  , 
and 

w  =  COS0  . 
In  order  to  find  scattered  fields,  the  current  must  be  known.  The  method  of 
moments  will  be  used  to  find  the  current.  In  this  approach,  the  unknown  quantity  is 
represented  by  a  series  with  unknown  expansion  coefficients.  The  boundary  conditions  are 
imposed  to  obtain  an  equation  for  the  current.  After  the  current  is  found  the  scattered  fields 
can  be  obtained  from  equations  (10)  and  (11).  This  is  referred  to  as  the  direct  problem,  i.e., 
finding  the  fields  from  the  currents.  Determining  the  currents  from  the  far  field  is  the 
inverse  or  synthesis  problem. 


D.         RCS  SYNTHESIS  APPROACH 

RCS  synthesis  is  a  part  of  a  more  general  problem  of  inverse  scattering.  In  an 
inverse  scattering  problem,  the  scattering  field  is  known  and  the  objective  is  to  extract 
current  and  resistivity  information  about  the  target  from  the  scattered  field.  There  are  two 
types  of  inverse  scattering  problems:  general  and  restricted.  In  the  general  case  nothing  is 
known  about  the  target  whereas  in  the  restricted  case  some  characteristics  of  the  target  are 
known.  The  radar  cross  section  synthesis  problem  falls  into  the  latter  case.  The  target  shape 
is  known,  but  the  surface  materials  and  current  are  to  be  determined.  Synthesis  problems 
can  be  solved  in  the  time  or  frequency  domains.  This  thesis  deals  with  synthesis  problems 
in  the  frequency  domain. 

At  an  observation  point  P,  the  scattered  electric  field  will  be  [3] 


Up)  =  -M  j[ 

where 


JsG  +  ^{h'xJs)xV'G-^-4^V'G\is'        (15) 


e~jkR 
G  =  —— —  and  R  = 

4nR 


r  —r 


Primed  quantities  are  associated  with  the  source  point;  unprimed  quantities  refer  to  the 
observation  point.  In  this  equation,  ES(P)  is  known,  but  the  current  and  surface 
impedance  are  both  unknown.  This  equation  can  be  solved  approximately  or  rigorously. 
In  the  first  case,  the  current  can  be  approximated  and  the  equation  can  be  solved  for 
impedance.  Equation  (15)  is  sufficient  to  completely  describe  the  problem.  In  the  second 
case,  both  current  and  impedance  must  be  determined.  This  requires  additional 
information  provided  by  the  following  equation:  [3] 
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^^  -  Z,(rW)  x  fll(F)  =  Z,(r)n(r)  x  [v  x  \[j fids' 
+  /ft  }|Z,(«'  x  Js)Gds'  +  j-  j J  V  ■  \Zsn'  x  J^Gais' 


+ 


-I  tan 


(16) 

/ft  JP.G*-  +  ^  J|(V  ■  £)VGcfc'  -  V  x  JJz,(»'  x  J.Jaft- 

Equations  (15)  and  (16)  can  be  solved  simultaneously  for  the  current  and  impedance. 
E.         METHOD  OF  MOMENTS 

The  method  of  moments  is  a  frequency  domain  technique.  It  is  used  for  solving 
complex  integral  equations  by  reducing  them  to  a  system  of  linear  equations.  The  MM 
technique  requires  that  the  entire  structure  to  be  modeled  be  broken  down  into  subdomains. 
The  surface  is  subdivided  to  a  number  of  patches  with  dimensions  small  compared  to  the 
wavelength.  Moment  methods  employ  a  technique  known  as  the  "method  of  weighted 
residuals."  Actually,  the  terms  method-of-moments  and  method-of-weighted-residuals  are 
synonymous  [5]. 

The  MM  technique  is  popular  because  it  can  be  applied  to  arbitrary  bodies. 
Furthermore,  integral  equations  are  reduced  to  a  set  of  linear  equations,  which  are 
conveniently  solved  by  using  matrix  algebra.  The  size  of  the  matrix  is  directly  related  to 
the  electrical  size  of  the  body.  Large  bodies  result  in  large  matrices,  which  affects 
computation  time,  and  limits  the  size  of  the  targets  that  can  be  handled  in  practice. 

The  first  step  in  MM  for  the  direct  problem  is  to  represent  the  current  by  a  series 
with  unknown  expansion  coefficients.  For  a  rectangular  resistive  plate  in  the  xy  plane,  one 
possible  representation  for  the  current  is: 

wy)=S?S+iS).  a?) 


ii 


The  vectors  J*  and  J£  are  the  expansion  or  the  basis  functions,  which  are  complex  in 
general.  Selection  of  the  basis  functions  is  very  important.  They  should  be 
mathematically  convenient  -  in  other  words,  easy  to  integrate  and  differentiate.  They 
should  also  be  consistent  with  the  behavior  of  the  current  to  provide  fast  convergence 
[3].  The  expansion  scheme  for  the  currents  can  be  seen  in  Figure  3.  Two  basis  functions 
span  each  rectangular  subsection:  one  x-directed  and  a  second  y-directed.  These  two 
orthogonal  basis  functions  are  sufficient  to  uniquely  define  all  possible  current  vectors  on 
the  patch. 


Figure  3.  Three-dimensional  View  of  Pulse  Basis  functions  for  the  Currents  Jx  and  J 
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III.  FORMULATION  AND  SOLUTION 


A.         RCS  OF  A  RESISTIVE  SHEET  USING  THE  METHOD  OF  MOMENTS 

In  this  section,  the  scattered  field  and  radar  cross  section  of  a  resistive  sheet  is 
calculated  for  the  bistatic  case.  This  direct  problem  solution  is  used  to  generate  and 
validate  the  synthesis  procedure.  A  constant  resistivity  is  assumed  on  each  patch,  and  the 
method  of  moments  is  used  to  find  the  current  distribution  on  the  plate.  Rectangular  patches 
and  nonoverlapping  pulse  basis  functions  are  used  for  calculations.  After  the  synthesis 
procedure,  the  resulting  resistivity  is  compared  with  the  original  resistivity  to  verify  the 
synthesis  equations. 

The  resistive  sheet  lies  in  the  xy  plane.  The  incident  wave  is  assumed  to  be  a  planar 
wave,  and  it  contains  both  parallel  and  perpendicular  components  [3], 


where, 


and, 


E'  =  (e^  +W*Yfk~r  =  (xE'x  +yEiyy*i 


-jk-r  =  jk[x s'mBj  coscj),  +  ^81110,  sin^  +  ZCOS0,] 


(18) 


(19) 


(20) 


El  =  Eq  cos6z-  cos(j)y  -  E'^  sin<t>z- 
E    =  Eq  cosG  •  sin<{>  •  +  E1^  cos<|>  • 

The  resistive  sheet  integral  equation  is  solved  to  find  the  current  distribution  on  the  surface 
[3], 


m 


=  RSJS  + 


tan 


jh]Q[(jsGdS'+^ 


V»jjJsGds' 


(21) 


tan 
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(22) 


(23) 


The  unknown  current  is  expanded  into  a  series  using  equation  (17),  which  can  be  written  in 
the  following  manner: 

Ay  Ax 

n  s  n 

The  quantities  Ax  and  Ay  are  defined  in  Figure  4,  and  pn(x\y)  is  the  pulse  function 

.        .      fl     ,  if  x'  and  y  on  patch  n. 
P"(x''y)  =  (o  .else 

The  method  of  moments  is  initiated  by  inserting  equation  (22)  into  (21).  At  this 
point  a  set  of  testing  or  weighting  functions  are  defined.  The  weighting  functions  describe 
how  the  boundary  conditions  are  enforced  [6].  These  functions  are  generally  chosen  to  be 

the  complex  conjugate  of  the  basis  functions  (  W  =  J  ),  which  is  known  as  "Galerkin's 
Method"  [1].  Weighting  functions  for  this  case  are  [7], 

W  =  Wxx  +  Wyy  (24) 

where, 

Wx=xL      7-  (25) 

m      y 


Wy=yiP">^ 


m 


Ax 


(26) 


Both  sides  of  the  equation  (21)  are  multiplied  with  these  weighting  functions  and  integrated. 
In  terms  of  expansion  and  weighting  functions,  the  integral  equation  (equation  21)  takes  the 
form 


jJ£,  •  Wds  =  l\RsJs  •  Wds  +  jkr\0  jjds  jjds' 


js.w-—(v.wp*.js) 


(27) 
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Equation  (27)  can  be  separated  into  two  equations,  one  associated  with  the  x 
components  and  the  second  with  the  y  components.  The  derivatives  of  the  pulse  functions 
are  approximated  by  5  functions  as  shown  in  Figure  4.  The  x  component  of  equation  (27) 
becomes, 


s(x-K,y) 


s{x,y-y~n) 


Figure  4.  Pulse  Function  Derivatives  for  a  Patch  [7]. 
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Ay         n  *[£  n    Ay  j    Ay 


s         v  n 


+jla\0\\ds\\ds<  \W 


Pn(x',y')  I  PmO'jO         1 


5  5 


Ay 


Ay  jfc' 


iM*-*-*) 


(28) 


-d(x-x+,yj}<  ^  ^L^x-X;,^_8^_x+ ,^  +  ^5^-^ J 


-s(x,7-y^)j  j  kr(x,x',>;,y) 


A  similar  equation  can  be  obtained  by  testing  with  the  y-directed  weighting  functions, 


5  V    « 


.     ...Ol 


+An.JWJ*   £^ 


pm(x,y) 1_ 

Ax  k2 


i;{i*>y-y~) 


(29) 


-B{x,y-yt)}  2  ^(5(x'-x„-,y)-8(x'-xM+,y))  +  £(5(x',y-^) 

-o(x\y-7:))]|lJG(x,x',7,y). 
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Now  define  the  following  quantities, 


s~iXX     _ 


Mi 

A/ 


7*10 


j J dxdy  \\dx 'dy'G(x,x',y,y') 


*m  °rt 


kAy' 


torn    ton 


\dy  \dy'\G{x-m,x-,y,y')-G[x-m,x+n,y,y') 


-Gix^,xn,y,y,]j  +  G^,x^y^y'j\ 


(30) 


(31) 


Qxy  _  -mp 

kAyAx 


jdy  ldx'\G\x~,x',y,yn)-G\xm,x',y,y+)  (32) 


Aym     Axn 


>R 


Mm        J 


Ax 
Rsm—    ,dm  =  n 

0         ,  if  m*  n 


(33) 


V;  =  E'x  Ax  sine 


'kAxuA  z^fkAyvAJk(XnUi+ynVi) 


G£=^jjdxdy\jdx'dy'G(x,x',y,y') 


Ax2 


sm  sn 


(34) 
(35) 


QZ=Tpr   fa  ^'\G(x,x',ym,yn)j-G(x,x',ym,y^ 


kAx' 


tym  4yw 


-G(x,x\y+,y„}  +  G{x,x\y+tytyi 


(36) 


17 


Axw     Ay„ 


(37) 


0    m 


Ay 


0 


if  m&  n 


V*  =  £'Aj>sinc 


kAxUi)  z^(k^i\jk(xnui+ynVi) 

K     2    J        \     2    ) 


(38) 


(39) 


Using  equations  (30)-(39),  equation  (28)  and  equation  (29)  can  be  cast  in  matrix  form  : 

f 


Vx 
Vy 


r 
r 


+ 


r 
r 


(40) 


gxx+qxx  Qxy 

For  the  synthesis  problem,  the  number  of  unknowns  in  this  equation  is  4  N  where  N  is  the 
number  of  expansion  terms.  There  are  2N     current  coefficients  and  2N     resistivity 

coefficients  in  general.  If  the  material  is  isotropic  then  f*  =  fv  and  the  number  of  unknown 
resistivity  coefficients  is  just  N .  For  the  direct  problem  the  number  of  unknowns  is  only 
2N. 

The  integrals  appearing  in  the  above  equations  must  be  evaluated  numerically.  To 
reduce  the  computer  run  time,  the  following  approximation  can  be  applied  [3], 

A 


xm  + 


jf(x,x')dxzAf(xm,x') 


(41) 


xm 


With  this  approximation,  the  impedance  matrix  elements  in  equations  (31),  (32),  (36)  and 
(37)  are 


Qmn  =  ^^{G(xm>xn>ym>yn)-G(xm>xZ>ym>yn) 

-G{x+,x~,ym,yn)  +  G\x+,x];fym,y„) 


(42) 
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QZ  =^^{G(xm>xn>ym>yn)-G(xm>xri>yn>yZ) 

-G\x+,xn,ym,y-)  +  G\x+,xn,ym,y+y> 

Qmn=-^\G[xm,x~,yl,ynYG{xm,x^,yl,yn^ 
-G\xm,x+,y-,yn)  +  G{xm,x+,y+,ynyi 

Q^=^^\G[xm,xn,y~,y~^-G^xm,xn,y~,y^ 
-G\xm,x„,y+,y-]  +  G\xm,x„,y+,y+ 


(43) 


(44) 


(45) 


Note  that  in  all  of  these  equations,  the  Green's  functions  in  the  integrands  can 
become  singular  for  the  self  term  because  both  the  source  and  the  test  points  lie  on  the  same 
segment.  This  problem  can  be  solved  by  moving  the  test  point  off  the  patch  axis  [3].  For 
example,  when  integrating  in  x  or  x ', 


R  =  ^(x-x')2+{y-y'f+a 


2 

y 


(46) 


and  when  integrating  in  v  orjy', 


R  =  j(X-X'f+(y-y>)2+ax: 
In  this  case,  an  equivalent  radius  of  the  patch  is  used  : 


(47) 


ax  =  0.225  Ax         ay  =  0.225  Ay  (48) 

Equation  (40)  provides  a  means  of  solving  the  direct  problem  by  plugging  the 
resulting  current  coefficients  into  the  far-field  radiation  integrals.  The  current  coefficients 
can  be  obtained  if  the  incident  field  and  resistivity  are  known.  The  scattered  field  at  an 
observation  point  (Qs  ,§s)  is  computed  by  integrating  the  current.  The  current  coefficients 
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that  occur  in  the  scattering  equations  are  the  same  for  all  observation  angles,  because  the 
same  incident  field  determines  these  coefficients.  After  integrating,  the  scattered  electric 
field  has  the  following  form: 


p   _  -jfa\p  c-jkr 


E   = 


+ 


(    N 


(     N 


.     ( kAxuA   .     fkAyv^ 
sine 


J  kg 


f k/Sxus 
\     2 


2    )        \     2    ) 
\ 


sine 


( kAyv 


\     2 


s_  \e-JkS 


J 


(49) 


xEsx+yEsy 


where: 


g  =  x„us+ynvs 
=  xn  sin05  cos<{>5  +yn  sin05  sin^ 

Finally,  the  complete  expression  for  the  0-  and  <j>-  components  of  the  scattered  electric  field 
is: 


Eq  -  Esx  cos05  coscj)  5  +  Esy  cos05  shuj^ 

Esif=-Esx$m§s+Esycos§s 

From  the  definition  of  the  radar  cross  section  for  a  three-dimensional  target  [1] 


(50) 
(51) 


cj=  lim 

r— >oo 


4nr 


m2 


(52) 


and  the  corresponding  radar  cross  section  for  the  0-  and  (j>-  components  of  the  scattered 
electric  field  for  a  unit  plane  wave  are 

k\20r       _,2 


aee  "      471 
^"      4n 


V 


(53) 
(54) 
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The  first  subscript  on  a  denotes  the  polarization  of  the  scattered  field,  whereas  the  second 
one  denotes  the  polarization  of  the  incident  wave. 

B.         SYNTHESIS  FORM  OF  THE  SCATTERING  EQUATIONS 

The  resistivity  of  the  rectangular  plate  can  be  determined  rigorously  by  expanding 
both  current  and  resistivity  into  series  and  solving  for  the  expansion  coefficients.  For  this 
purpose,  the  scattered  electric  field  is  found  using  the  formulation  in  the  previous  section. 
Because  the  sine  functions  in  equations  (34)  and  (39)  are  periodic  in  direction  cosine  space 
(DCS),  it  is  convenient  to  work  in  DCS  (Figure  5).  There  are  two  approaches  sampling  the 
RCS  in  DCS.  The  first  is  to  choose  the  observation  points  entirely  inside  of  the  unit  circle. 
The  second  is  to  include  observation  points  outside  of  the  unit  circle.  The  first  approach 
makes  physical  sense  because  all  observation  directions  are  points  in  real  space.  Directions 
outside  of  the  unit  circle  correspond  to  complex  angles.  However,  the  first  method  is 
restrictive  in  the  number  of  observation  points,  because  they  must  be  equal  to  the  number  of 
subdomains  on  the  plate  to  obtain  a  square  matrix  that  can  be  solved  by  inversion.  The 
computer  program  which  calculates  and  plots  RCS  in  DCS  is  presented  in  the  Appendix. 

In  synthesizing  the  plate  resistivity  the  far-field  pattern  is  given.  It  is  known  that  a 
resistive  sheet  does  not  support  magnetic  currents,  and  the  far-field  is  related  to  the  electric 
currents  on  the  sheet  by 

Es(P)  =  -jkr\0jjJsGds'  (55) 

5 

P  denotes  an  observation  point  in  the  direction  (r,6,<j>)  or  {x,y,z)  .  On  the  surface  of  the 
sheet,  equation  (21)  still  holds,  and  in  this  equation  resistivity  and  current  are  not 
independent  of  each  other.  These  two  quantities  are  related  at  every  point  on  the  surface  of 
the  sheet  by  the  tangential  component  of  the  electric  field  as  indicated  in  equation  (4). 

Equations  (21)  and  (55)  can  be  solved  simultaneously  for  the  current  and  resistivity. 
For  the  current  expansion  equation  (22)  is  used,  and  for  the  resistivity, 
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RS  =Y.RnPn{* \y') 


(56) 


u 


Unit  Sphere 


►  v 


+  V  =  1 


Figure  5.  Direction  Cosine  Space  (DCS) . 
The  resistivity  is  assumed  to  be  constant  on  each  patch.  Equation  (22)  is  plugged  into  the 
equation  (55),  and  after  integration, 


Es{Pm)  =  -M 

where  the  last  factor  in  the  equation  is 


-jkr    N  f  rx  ry 

47tr  ^TjvAy        Ax 


Fn(um,vm)  (57) 


Fn{um>vm)  =  ejk(XnUm+ynVm)  AxAysmci 


(  kAxu 


f  1r\yn>     ^ 


m 


V     2 


sine 


&Ayv 

k   2   ; 


(58) 
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Equation  (57)  can  be  solved  for  the  unknown  current  coefficients  given  a  set  of  m  = 
1,2,.... ,N  scattered   field  values.  In  matrix  form,  equation  (57)  is 


Ej=— I*        Eys=—ly 
Ay  Ax 

These  two  equations  can  be  combined  and  put  into  the  following  form 


(59) 


_            _ 

Ay 

E* 

lx 

_E?_ 

F/ 

y 

(60) 


In  equation  (60),  the  only  unknowns  are  the  I*  and  V  coefficients,  and  they  can  be  found 
by  simple  matrix  algebra.  After  these  coefficients  are  found,  they  are  plugged  into  the 
resistive  sheet  integral  equation  (equation  (21))  to  find  the  resistivity  coefficients  in 
equation  (56). 


Vy 
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r 
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Zxx 

Zyx 


zxy 

Zyy 
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(61) 


There  are  4N  rows  in  equation  (61).  In  an  isotropic  resistivity  case,  there  are  3  N 
unknowns  because  the  terms  p  x  and  p  y  in  equation  (61)  are  same.  Rearranging  equation 
(61) 


fr" 
iy 

= 

vy 
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Z" 

Zyx 


zxy 

Zyy 


r 
r 


(62) 


The  matrix  operations  on  the  right  hand  side  yield  a  column  vector.  The  elements  of  the 
column  vector  are  fxmVm  and  ^  yJym  and  since  /*  and  Iym  are  known,  the  resistivity  can 
be  found  by  element-wise  division. 
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IV.       COMPUTER  IMPLEMENTATION  AND  RESULTS 

A.  INTRODUCTION 

In  this  chapter,  the  computer  coding  of  the  RCS  synthesis  formulas  is  described  and 
numerical  results  are  presented  for  a  flat  resistive  sheet.  The  results  of  three  investigations 
are  included : 

1.  Verification  of  the  MM  solution  for  rectangular  patches  and  pulse  basis 
functions. 

2.  Verification  of  the  synthesis  equations  derived  in  Chapter  III  and  a  parametric 
study  of  convergence. 

3.  Sampling  requirements. 

B.  VALIDATION  OF  THE  RECTANGULAR  PATCH  CODE 

The  rectangular  patch  code  based  on  pulse  expansion  functions  developed  in 
Chapter  III  was  programmed  in  MATLAB.  Computed  results  were  compared  to  a  MM 
code  named  PATCH  [9].  Typical  results  are  shown  in  Figures  6  through  10.  The  two  codes 
give  essentially  the  same  RCS  as  expected.  The  agreement  is  within  2  dB  down  to  the  -40 
dB  level  when  the  two  codes  are  converged. 

Validation  of  the  rectangular  patch  code  is  important  because  many  of  the  same 
subroutines  and  functions  are  used  in  the  synthesis  program.  Furthermore,  the  output  of  this 
program  is  used  as  the  input  for  initializing  the  synthesis  process.  Figures  1 1  and  12  show 
two  dimensional  plots  of  the  RCS.  The  corresponding  electric  field  values  are  used  to  fill 
the  vector  V  in  Equation  (40). 

C.  VERIFICATION  OF  THE  SYNTHESIS  EQUATIONS 

The  synthesis  equations  and  computer  codes  are  validated  by  the  following  process  : 

1 .    Compute  the  scattered  field  for  a  plate  with  known  surface  resistivity. 
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2.  Use  the  scattered  field  from  step  1  in  the  synthesis  procedure. 

3.  Compare  the  synthesized  resistivity  in  step  2  to  the  original  resistivity  used  in 
step  1.  In  principle  the  difference  should  be  zero.  However,  because  of 
numerical  and  sampling  errors,  the  residual  will  not  be  exactly  zero,  but  some 
small  (and  hopefully  negligible)  value. 

The  difference  between  the  original  and  synthesized  resistivity  ARS  is  defined  by 

&RS  =  Rsy~  Rso 
where 

Rso  =  Original  resistivity  used  to  compute  the  electric  field  values. 
Rsy  =  Synthesized  resistivity. 

For  Rso=  377,  ARS  <  10"12  as  shown  in  Figure  13.  For  Rs=  10,  ARS  is  as  large  as  0.06 
as  illustrated  in  Figure  14.  Similar  results  are  shown  in  Figures  15  through  19.  The  lower 
plate  resistivity  supports  surface  waves  which  increase  the  size  of  the  off-diagonal  elements 
in  Z  relative  to  the  size  of  the  diagonal  elements.  This  leads  to  increased  roundoff  error.  The 
number  of  field  (RCS)  samples  Nu  and  jVv  are  equal,  and  they  include  the  points  outside 
of  the  visible  region,  which  is  the  second  approach  as  described  in  Chapter  III-B.  The 
number  of  patches  are  Nx  =  Nu  and  Ny  =  Nv.  All  of  the  visible  region  field  points  were 

in  the  backscatter  (rear)  hemisphere.  Including  the  points  outside  the  visible  region  affects 
the  accuracy  of  the  computation  and  causes  some  poorly  conditioned  matrices.  However, 
the  error  resulting  from  this  approach  is  not  significant.  In  all  cases  the  residual  resistivity 
was  less  than  0.6  percent. 

D.         SAMPLING  REQUIREMENTS 

Computer  memory  and  execution  time  are  very  important  limitations  in  the 
scattering  and  synthesis  computations.  Reducing  computer  run  time  depends  on  reducing 
the  size  of  the  matrices  (especially  impedance  matrix),  and  consequently  minimizing  the 
number  of  patches.  Thus  a  guideline  for  the  minimum  number  of  subdomains  which  gives 
computationally  accurate  results  is  necessary  to  minimize  the  computer  run  time. 
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For  the  synthesis  calculation  the  continuous  scattered  field  is  sampled  just  as  a 
continuous  time  function  is  sampled  in  a  digital  communication  system.  Therefore  the  same 
sampling  requirements  apply.  If  the  frequency  content  of  the  scattered  field  is  plotted,  the 
minimum  sampling  rate  must  be  at  least  twice  the  highest  frequency  [3].  To  obtain  the 
frequency  content  a  two-dimensional  discrete  Fourier  transform  (DFT)  is  performed  on  the 
scattered  electric  field  values,  which  are  computed  in  sufficiently  large  number.  The 
computation  is  performed  in  DCS  because  most  RCS  patterns  are  highly  periodic  in  DCS 
and  thus  the  number  of  sampling  points  is  smaller.  If  Nu  and  Nv  are  the  number  of 
sampling  points  in  the  u  and  v  coordinates,  then  the  sampling  intervals  are  2/Nu  and  2/Nv. 
This  represents  n  radians  in  visible  space,  and  hence  the  scaling  of  the  u  and  v  axes  upon 
returning  from  the  FFT  is  2/Nu  Au  and  2/Nv  Av ,  respectively. 

In  Figures  20  through  23,  frequency  information  for  Yk  x  Yk  and  3k  x  3A.  plates  is 
plotted.  The  FFT  is  performed  on  the  center  row  of  the  two-dimensional  scattered  electric 
field  matrix  to  illustrate  the  frequency  content  in  one  coordinate.  The  magnitude  of  this 
transform  is  plotted  with  respect  to  the  number  of  cycles  in  DCS  in  Figures  20  and  22.  The 
result  of  a  two-dimensional  FFT  performed  on  the  whole  E-field  matrix  is  plotted  with 
respect  to  the  number  of  cycles  in  DCS  in  Figures  21  and  23.  The  maximum  frequency  for 
Yk  x  Yk  resistive  plate  is  1  for  the  bistatic  case.  Since  the  sampling  frequency  must  be  twice 
the  maximum  frequency  component  of  the  signal,  the  sampling  frequency  is  2,  which 
corresponds  to  0.5A.  x  0.5A,  patch  size.  Applying  this  procedure  to  a  3k  x  3A,  plate  yielded 
the  same  maximum  patch  size.  The  required  sampling  frequency  approaches  zero  as  the 
plate  size  gets  smaller,  as  expected.  Thus  in  the  case  of  bistatic  RCS  synthesis  a  half- 
wavelength  subdomain  size  corresponds  to  the  Nyquist  sampling  rate.  Note  that  this 
guideline  applies  only  to  cases  where  the  RCS  pattern  is  highly  periodic  (i.e.,  negligible 
multiple  reflections,  surface  waves,  diffraction,  etc.).  When  these  other  scattering 
mechanisms  are  present  the  scattering  pattern  becomes  more  complicated  and  thus  contains 
more  high-frequency    components.  Also,  it  is  important  to  note  that  satisfying  the  field 
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sampling  rate  does  not  insure  that  the  series  expansion  for  the  current  on  the  plate  is 
converged. 
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Figure  6.     Comparison  of  the  bistatic  RCS  computed  with  PATCH  vs  that  with 

rectangular  patches  (Rs  =  377,  Nx=Ny  =  8,  Lx=  Ly=  lA  and  #/=  <j)f=  0°). 
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Figure  7.     Comparison  of  the  bistatic  RCS  computed  with  PATCH  vs  that  with  rectangular 
patches  (Rs=  377,  Nx  =  Ny=  14,  Lx  =  Ly  =  U and  9=  ft=  0°). 
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.80      -60      -40      -20         0        20       40        60        80 
<?,  degrees  (<j>  =  0°) 


Figure  8.     Comparison  of  the  bistatic  RCS  computed  with  PATCH  vs  that  with  rectangular 
patches  (Rs=  10,  Nx  =  Ny=  20,  Lx  =  Ly=  U  and  &=  ft=  0°). 
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<?,  degrees  (<j>  =  0°) 

Figure  9.     Comparison  of  the  bistatic  RCS  computed  with  PATCH  vs  that  with  rectangular 
patches  (Rs  =  377,  JVX  =  JVy  =  20,  Lx  =  Ly  =  U  and  6X  =  30°,  $=  0°). 
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-80      -60      -40      -20        0        20       40       60       80 
<?,  degrees  (<j>  =  0°) 

Figure  1 0.  Comparison  of  the  bistatic  RCS  computed  with  PATCH  vs  that  with  rectangular 
patches  (Rs=  377, Nx  =  Ny=  20,  Zx=  L=  21  and  &=  fr=  0°). 
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v  =  sin  <?sin  <f> 
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&  =  sin  tfcos^ 


Figure  1 1 .  Two  dimensional  plot  of  the  bistatic  RCS  of  a  rectangular  resistive  sheet  (Rs: 
377,  Nu  =  Nv=  23,  (Nx  =  Nv ,  Ny  =  NJ,  L=  L=  U  and  fy  =  fr  0°). 
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Figure  12.  Two  dimensional  plot  of  the  bistatic  RCS  of  a  rectangular  resistive  sheet  (Rs 
377,  Lx=  Ly=  U  and  <9/=  30°,  $=  0°). 
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Figure  13.  Residual  resistivity  for  Rso  =  377,  N^Ny^  (N^NU,  Ny=Nv),  Lx=  Ly=  1A  and 
normal  incidence. 
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Figure  14.  Residual  resistivity  for  Rso  =  10,  Nu = JVV=  9  (A^=  iVu,  Mf=  NvX  ^x  =  £y=  1  ^ 
and  normal  incidence. 


36 


-1    -1 


Figure  1 5 .  Residual  resistivity  for  R so  =  377,  NU  =  NV=  16  (Nx=  Nu,  Ny=  Nv),  Lx  =  Ly=  1 A 
and  normal  incidence. 
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Figure  16.  Residual  resistivity  for  Rso  =  377,  Nu  =  Nv  =  1 6  (Nx=  Nu,  Ny=  Nv),  Lx  =  Ly=  1 A 
and  $i=  30°,  0/=O°. 
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Figure  1 7.  Residual  resistivity  for  Rso  =  377,  Nu = Nv  =  1 6  (A^=  Nw,  Ay=  Nv),  Lx  =  Ly=  2/1 
and  normal  incidence. 
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Figure  1 8.  Residual  resistivity  for  Rso  =  377,  NU=NV=  23  (Nx=  Nu,  Ny=  Nv),  Lx  =  Ly=  1 X 
and  normal  incidence. 
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Figure  19.  Residual  resistivity  for  Rso=  377,  NU=NV=  23  (Nx=  Nu,  Ny=  Nv),  Lx=Ly=\X 
and  0j  =  30°,  #=0°. 
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Figure  20.   Frequency  analysis  of  one  row  (or  u  coordinate)  of  the  scattered  E-field.  Maximum 
frequency  is  1  for  7^=  377,  Lx=Ly-  \X,  256  field  points  and  normal  incidence. 
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Figure  21 .  Frequency  analysis  in  two-dimensions  for  the  scattered  E-field.  Maximum 
frequency  is  1  for  Rso=  377,  Lx=  Ly=  1A,  256  field  points  and  normal 
incidence. 
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Figure  22.   Frequency  analysis  of  one  row  (or  one  column)  of  the  scattered  E-field.  Maximum 
frequency  is  3  for  Rs0=  377,  Lx=  Ly=  3 A,  900  field  points  and  normal  incidence. 
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Figure  23.  Frequency  analysis  in  two-dimension  for  the  scattered  E-field.  Maximum 
frequency  is  3  for  Rs0  =  377,  Lx=  Ly=  3/1,  900  field  points  and  normal 
incidence. 
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V.  CONCLUSIONS 

A  rigorous  historic  RCS  synthesis  method  for  flat  rectangular  sheets  has  been 
presented.  For  convenience,  the  research  was  limited  to  resistive  surfaces  and  parallel 
polarization.  The  scattering  and  synthesis  equations  are  based  on  the  method  of  moments 
(MM)  technique.  A  rigorous  method  requires  that  both  the  current  and  resistivity  be 
expanded  into  series,  and  the  unknown  coefficients  determined  simultaneously. 

The  scattered  field  for  a  resistive  sheet  was  first  calculated  for  a  known  resistivity. 
The  historic  RCS  of  the  resistive  sheet  was  used  to  validate  the  rectangular  patch  solution 
by  comparing  it  with  the  output  of  PATCH  for  different  cases.  This  direct  problem  solution 
was  used  to  generate  and  validate  the  synthesis  procedure.  The  unknown  current  and 
resistivity  coefficients  were  found  by  solving  two  coupled  integral  equations  (equations  (15) 
and  (16)  ).  The  reconstructed  and  original  resistivities  were  compared  to  verify  the 
synthesis  equations.  A  parametric  study  of  the  solution  convergence  was  performed.  In  all 
cases,  the  reconstructed  resistivity  obtained  from  the  synthesis  equations  was  in  agreement 
with  the  original  resistivity. 

The  frequency  content  of  the  scattered  electric  field  was  also  investigated  in  order  to 
determine  the  maximum  allowable  patch  size.  The  number  of  far-field  points  must  be 
equal  to  the  number  of  patches  on  the  plate  to  obtain  a  unique  solution.  If  the  field  pattern  is 
highly  periodic  in  DCS,  as  in  the  case  of  the  physical  optics  solution,  the  Nyquist  sampling 
theorem  yields  approximately  0.5A,  patch  size  on  the  plate  for  bistatic  RCS.  However,  a 
MM  solution  with  0.5A.  patches  is  usually  not  converged;  the  conventional  guideline  is 
0.1A.  subdomains.  However,  in  many  cases  of  practical  interest  where  the  patterns  to  be 
synthesized  are  primarily  determined  by  specular  scattering  (i.e.,  negligible  surface  waves 
and  multiple  reflection  and  diffraction)  the  rninimum  patch  size  can  possibly  be  larger  than 
the  usual  0.1  A, . 

The  rigorous  solution  presented  in  this  thesis  is  limited  to  flat  rectangular  resistive 
sheets.  This  allowed  the  use  of  rectangular  subdomains  with  pulse  basis  functions.  Pulse 
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basis  functions  are  convenient  for  computer  programming  and  result  in  matrix  synthesis 
equations  that  are  easy  to  solve.  However,  the  more  important  scattering  problems  involve 
doubly  curved  surfaces.  For  these  targets  triangular  subdomains  must  be  used  to  accurately 
represent  the  surface.  The  overlapping  basis  functions  of  the  type  used  in  PATCH  will 
result  in  much  more  complicated  (and  possibly  nonlinear)  sets  of  equations.  It  is 
recommended  that  this  be  the  direction  of  future  research  in  this  area. 
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APPENDIX  COMPUTER  CODES 

In  this  appendix  computer  codes  are  presented.  The  code  plate,  m  calculates  the 
RCS  in  O-direction  for  the  bistatic  case  and  T1VL,  polarization.  The  code  platedcs.m 
calculates  the  RCS  and  plots  it  in  direction  cosine  space.  The  data  from  this  code  is  saved 
to  be  used  in  the  synthesis  program.  The  code  synt.m  uses  the  stored  data  to  reconstruct  the 
surface  resistivity.  The  codes  dcsin.m,  center. m,  green.m,  greenl.m,  green2.m,  green3.m, 
and  green4.m  are  the  functions  which  are  used  in  the  synthesis  and  scattering  programs. 
The  code  sample.m  performs  the  FFT  analysis  in  the  u  and  v  directions  of  the  scattered 
field.  Detailed  explanations  are  presented  in  the  comments  in  the  codes. 
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%%%%%%%%%%%%0/o%%0/o0/o%0/o0/o0/o0/o0/o0/o0/o0/o0/o%%%%0/o0/o0/o%%%%%%% 

%                                           plate.m  % 

%  Resistive  Square  Plate  Bistatic  RCS  (in  teta  direction)  % 

%  with  Method  of  Moments.  % 

%  % 

%  Ugurcan  Samli,  July  1996  % 

%  % 

%%%%%%%%%%%%%%%%%%%0/o0/o0/o0/o0/o0/o0/o0/o%0/o%%0/o0/o%%%%%%%% 

clear  all; 
format  compact; 

%%%%%%%%%%%%%%  Variables  used  in  the  program%%%%%%%%%%%%% 

%  modex  and  modey  are  the  number  of  patches  in  x-  and  y-directions 
%  respectively.  Eotet  and  Eophi  are  for  polarization.  Rs  is  the 
%  constant  surface  resistivity. 

Lx=l ;  %  input('Length  in  x-direction= '); 

Ly=  1 ;  %  input('Length  in  y-direction= '); 

modex=20;  modey=20; 
deltax=Lx/modex;  deltay=Ly/modey; 
mode=modex*modey;  eta=377;  lambda=l; 
k=2*pi;  t=pi/180;  Rreal=377;  Rimg=0; 
Eotet=l;  Eophi=0;  Emag=sqrt(EotetA2+EophiA2); 
incr=2;  sta=-90;  stp=90; 
lim=floor((stp-sta+ 1  e-5)/incr)+ 1 ; 
int=2;  phiang=0;  phi=phiang*t; 

Rsl=Rreal+j*Rimg;  Rs2=Rsl*eye(mode);  Rs=diag(Rs2)'; 

%%%%%%%%%  Computing  patch  center  points  %%%%%% 

in=0; 

for  s=l:modey 

for  ss=l:modex 

in=in+l; 

xn(in)=(-Lx/2+deltax/2+(ss- 1  )*deltax); 
yn(in)=(-Ly/2+deltay/2+(s- 1  )*deltay); 

end; 
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end; 

xm=xn;  ym=yn; 

funl=j*k*eta;  fiin2=j/k*eta; 

%%%%%%%%%%  Calculation  of  impedance  matrix  %%%%%%%%%%% 

%  Impedance  matrix  calculation  is  same  with  the  impedance  matrix 
%  calculation  in  platedcs.m,  green  functions  are  used  in  calculations. 

forn=l:mode 
form=l:mode 

Gxx=green(xm(m),xn(n),ym(m),yn(n),deltax,deltay,int); 
Gyy=green(ym(m),yn(n),xm(m),xn(n),deltay,deltax,int); 

xnl  =xn(n)+deltax/2; 
xn2=xn(n)-deltax/2; 
xm  1  =xm(m)+deltax/2; 
xm2=xm(m)-deltax/2 ; 
yn  1  =yn(n)+deltay/2; 
yn2=yn(n)-deltay/2 ; 
ym  1  =ym(m)+deltay/2 ; 
ym2=ym(m)-deltay/2; 

Qxx=(green  1  (xm2,xn2,ym(m),yn(n),deltax,deltay,int). . . 

-greenl  (xm2,xnl  ,ym(m),yn(n),deltax,deltay,int)... 

-greenl  (xml  ,xn2,ym(m),yn(n),deltax,deltay,int)... 

+green  1  (xm  1  ,xn  1  ,ym(m),yn(n),deltax,deltay  ,int)); 
Qxy=(green2(xm2,xn(n),ym(m),yn2,deltax,deltay,int). . . 

-green2(xm2,xn(n),ym(m),yn  1  ,deltax,deltay  ,int). . . 

-green2(xm  1  ,xn(n),ym(m),yn2,deltax,deltay  ,int). . . 

+green2(xm  1  ,xn(n),ym(m),yn  1  ,deltax,deltay  ,int)); 
Qyx=(green3(xm(m),xn2,ym2,yn(n),deltax,deltay,int)... 

-green3  (xm(m),xn2,ym  1  ,yn(n),deltax,deltay  ,int). .  . 

-green3  (xm(m),xn  1  ,ym2,yn(n),deltax,deltay  ,int). . . 

+green3(xm(m),xnl  ,yml  ,yn(n),deltax,deltay,int)); 
Qyy=(green4(xm(m),xn(n),ym2,yn2,deltax,deltay  ,int). .  . 

-green4(xm(m),xn(n),ym2,ynl,deltax,deltay,int)... 

-green4(xm(m),xn(n),yml,yn2,deltax,deltay,int)... 
+green4(xm(m),xn(n),yml  ,ynl  ,deltax,deltay,int)); 
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sum=fun  1  *  Gxx-fun2  *  Qxx; 
if  m=n,  sum=sum+Rs(m)*(deltax/deltay);  end 
Z(m,n)=sum; 

sum=-fun2  *  Qxy ; 
Z(m,n+mode)=sum; 

sum=-fun2  *  Qyx; 
Z(m+mode,n)=sum; 

sum=flinl  *Gyy-fun2*Qyy; 

if  m+mode==n+mode,  sum=sum+Rs(m)*(deltay/deltax);  end 
Z(m+mode,n+mode)=sum; 

end; 
end; 

ZI=pinv(Z); 

disp('Impedence  Matrix  Computed  and  Inverted'); 

%%%%%%%%%%%%%  Incidence  angle  %%%%%%%%%%% 

phii=0;  teti=0; 

sinphii=sin(phii*t);  sinteti=sin(teti*t); 

cosphii=cos(phii*t);  costeti=cos(teti*t); 

wi=sinphii*sinteti; 
uui=sinteti  *  cosphii; 

Ex=Eotet*  costeti  *  cosphii-Eophi  *  sinphii; 
Ey=Eotet*  costeti  *  sinphii+Eophi*  cosphii; 

for  n=l:mode; 
argx=k*deltax*uui/2; 
argy=k*deltay*wi/2; 
sincx=l;  sincy=l; 

if  abs(argx)>le-5,  sincx=sin(argx)/argx;  end 
if  abs(argy)>le-5,  sincy=sin(argy)/argy;  end 
expxy=exp(j*k*(xn(n)*uui+yn(n)*wi)); 
vx(n)=Ex*deltax*expxy*sincx*sincy; 
vy  (n)=Ey  *  deltay  *  expxy  *  sincx  *  sincy ; 
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end; 

%  Excitation  vector  V,  and  current  coefficients  I. 

V(l:mode)=vx(:); 
V(mode+l  :2*mode)=vy(:); 
I=ZI*V.'; 

%%%%%%%%  calculation  of  the  scattered  fields  %%%%%%%%% 

for  aa=l:lim; 

angle(aa)=(aa- 1  )*incr+sta; 

teta=angle(aa)*t; 

costet=cos(teta); 

sintet=sin(teta); 

cosphi=cos(phi); 

sinphi=sin(phi); 

w=sinphi*  sintet; 

uu=sintet*cosphi; 

forn=l:mode 

argx=k*deltax*uu/2; 

argy=k*deltay*w/2; 

sincx=l;  sincy=l; 

if  abs(argx)>le-5,  sincx=sin(argx)/argx;  end 

if  abs(argy)>le-5,  sincy=sin(argy)/argy;  end 

expxy=exp(j  *k*(xn(n)*uu+yn(n)*  w)); 

vx  1  (n)=deltax*  expxy  *  sincx*  sincy ; 

vy  1  (n)=deltay  *  expxy  *  sincx*  sincy ; 
end; 

Exs=vxl  *I(1  :mode); 
Eys=vy  1  *I(mode+l  :2*mode); 

%%%%  TM  and  TE  polarization  cases  are  together.%%%%% 

etet=-j  *k*  eta/2/k*  (Exs*  costet*  cosphi+Ey  s*  costet*  sinphi); 

ephi=-j*k*eta/2/k*(-Exs*sinphi+Eys*cosphi); 

etetrcs(aa)=etet; 

ephircs(aa)=ephi; 

end; 
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%%%%%%%  Calculation  of  the  RCS,  and  plots  %%%%%%% 

RCStet=10*loglO(abs(etetrcs).A2*2*k+le-10); 
RCSphi=10*loglO(abs(ephircs).A2*2*k+le-10); 

plot(angle,RCStet),grid; 

max  1  =max(max(RCStet)); 

db=(floor(max(maxl)/5)+l)*5; 

axis([-90,90,db-50,db]); 

xlabel('Angle  in  degrees'); 

ylabel('RCS  in  dBsm '); 

title("Bistatic  RCS  of  a  Resistive  Plate  With  MM'); 

%save  fl  Z  RCStet  angle  Rreal 
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%%%%%%%%%%%%%%%%%%%%%0/o0/o%0/o0/o%%0/o%%%%0/o0/o0/o0/o%%%% 

%  platedcs.m                                                                 % 

%  This  program  calculates  Resistive  Square  Plate  Bistatic  RCS  with                       % 

%  Moment  Method,  and  it  plots  this  RCS  in  direction  cosine  space.                         % 

%  % 

%  Ugurcan  Samli.                                                                                                  % 

%  % 

%%%%%%%%%%%%%%0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o%0/o0/o0/o%0/o0/o0/o%0/o0/o0/o%% 

clear  all; 
format  compact; 

%%%%%%%%%%%%  Variables  used  in  the  program  %%%%%%%%%%% 

Lx=l ;  %  input('Length  in  x-direction= '); 

Ly=l ;  %  input('Length  in  y-direction= '); 

T=  1 ;  %  If  the  number  of  field  points  equal  to  number  of 

%  patches,  T=l. 
%  modex  and  modey  are  the  number  of  segments  in  x  and  y  direction. 
%  Eotet  and  Eophi  determines  the  polarization  (Eotet=l  and  Eophi=0 
%  is  the  Teta  polarization.) 
eta=377;  lambda=l;  k=2*pi;  t=pi/180;  int=2; 
Rreal=377;  Rimg=0; 
Eotet=l;  Eophi=0; 

%  Finding  center  of  the  patches  in  Direction  Cosine  Space. 

%  dcsin.m  and  center.m  functions  are  called. 

modex=4;  modey=3; 

deltax=Lx/modex;  deltay=Ly/modey; 

mode=modex*modey; 

load  uvl  %  u  and  v  points  are  computed  seperately. 

%  Constant  surface  resistivity  Rs  matrix. 
Rsl=Rreal+j*Rimg;  Rs2=Rsl*eye(mode);  Rs=diag(Rs2)'; 

%%%%%%%%%%  Computing  Patch  center  points.  %%%%%%%%%%%%% 

in=0; 

for  s=l:modey 

for  ss=l:modex 

in=in+l; 

xn(in)=(-Lx/2+deltax/2+(ss- 1  )*deltax); 
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yn(in)=(-Ly/2+deltay/2+(s- 1  )*deltay); 
end; 
end; 

xm=xn;  ym=yn; 

funl=j*k*eta;  fun2=j/k*eta; 

%%%%%%%  Filling  Impedence  Matrix  %%%%%%%%%%%% 

%  Green  functions  are  called  in  the  impedance  matrix  calculation. 
forn=l:mode 
form=l:mode 

Gxx=green(xm(m),xn(n),ym(m),yn(n),deltax,deltay,int); 
Gyy=green(ym(m),yn(n),xm(m),xn(n),deltay,deltax,int); 

xn  1  =xn(n)+deltax/2 ; 
xn2=xn(n)-deltax/2; 
xm  1  =xm(m)+deltax/2 ; 
xm2=xm(m)-deltax/2; 
yn  1  =yn(n)+deltay/2; 
y  n2=y  n(n)-deltay  12 ; 
ym  1  =ym(m)+deltay/2 ; 
ym2=ym(m)-deltay/2 ; 

Qxx=(greenl(xm2,xn2,ym(m),yn(n),deltax,deltay,int)... 

-greenl  (xm2,xnl  ,ym(m),yn(n),deltax,deltay,int)... 

-green  1  (xm  1  ,xn2,ym(m),yn(n),deltax,deltay,int). .  . 

+greenl  (xml  ,xnl  ,ym(m),yn(n),deltax,deltay,int)); 
Qxy=(green2(xm2,xn(n),ym(m),yn2,deltax,deltay,int)... 

-green2(xm2,xn(n),ym(m),ynl,deltax,deltay,int)... 

-green2(xml,xn(n),ym(m),yn2,deltax,deltay,int)... 

+green2(xml  ,xn(n),ym(m),ynl  ,deltax,deltay,int)); 
Qyx=(green3  (xm(m),xn2,ym2,yn(n),deltax,deltay  ,int). . . 

-green3(xm(m),xn2,yml,yn(n),deltax,deltay,int)... 

-green3  (xm(m),xn  1  ,ym2,yn(n),deltax,deltay  ,int) . . . 

+green3  (xm(m),xn  1  ,ym  1  ,yn(n),deltax,deltay  ,int)); 
Qyy=(green4(xm(m),xn(n),ym2,yn2,deltax,deltay,int). . . 

-green4(xm(m),xn(n),ym2,ynl,deltax,deltay,int)... 

-green4(xm(m),xn(n),yml,yn2,deltax,deltay,int)... 

+green4(xm(m),xn(n),yml  ,ynl  ,deltax,deltay,int)); 
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sum=fun  1  *  Gxx-fiin2  *  Qxx; 
if  m==n,  sum=sum+Rs(m)*(deltax/deltay);  end 
Z(m,n)=sum; 

sum=-fun2*Qxy; 
Z(m,n+mode)=sum; 

sum=-fun2  *  Qyx; 
Z(m+mode,n)=sum; 

sum=fiinl  *Gyy-fun2*Qyy; 

if  m+mode=n+mode,  sum=sum+Rs(m)*(deltay/deltax);  end 
Z(m+mode,n+mode)=sum; 


end; 


end; 


ZI=inv(Z); 

disp('Impedence  Matrix  Computed  and  Inverted'); 

%%%%%%%  Incident  Angles  %%%%%%%%% 

phii=0;  teti=0; 

sinphii=sin(phii*t);  sinteti=sin(teti*t); 

cosphii=cos(phii*t);  costeti=cos(teti*t); 

wi=sinphii  *  sinteti ; 
uui=sinteti  *  cosphii ; 

Ex=Eotet*  costeti  *  cosphii-Eophi  *  sinphii; 
Ey=Eotet*  costeti  *  sinphii+Eophi  *  cosphii; 

forn=l:mode; 
argx=k*deltax*uui/2; 
argy=k*deltay*wi/2; 
sincx=l;  sincy=l; 

if  abs(argx)>le-5,  sincx=sin(argx)/argx;  end; 
if  abs(argy)>le-5,  sincy=sin(argy)/argy;  end; 
expxy=exp(j*k*(xn(n)*uui+yn(n)*wi)); 
vx(n)=Ex*deltax*  expxy  *  sincx*  sincy ; 
vy(n)=Ey  *  deltay  *expxy  *  sincx*  sincy; 
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end; 

%  Excitation  vector  and  current  coefficients,  V  and  I  respectively. 

V(l:mode)=vx(:); 

V(mode+l  :2*mode)=vy(:); 

I=ZI*V.*; 

%%%%%%%%%  calculation  of  the  field  points%%%%%%%% 

form=l:mode 

unitl  (m)=(sqrt(u(m)A2+v(m)A2)); 
phi(m)=(atan2(u(m)+l  e-7,v(m))); 
teta(m)=asin(unitl  (m)); 
costet=cos((teta(m))) ; 
sintet=sin((teta(m))) ; 
cosphi=cos(phi(m)); 
sinphi=sin(phi(m)); 

forn=l:mode 

argx=k*deltax*u(m)/2; 

argy=k*  deltay  *  v(m)/2 ; 

sincx=l;  sincy=l; 

if  abs(argx)>le-5,  sincx=sin(argx)/argx;  end 

if  abs(argy)>le-5,  sincy=sin(argy)/argy;  end 

expxy=exp(j*k*(xn(n)*u(m)+yn(n)*v(m))); 

vxl  (n)=deltax*expxy*sincx*sincy; 

vy  1  (n)=deltay  *expxy  *  sincx*  sincy ; 
end; 

%  x-  and  y-components  of  the  scattered  field. 
Exs=vxl*I(l:mode);  Exsl(m)=Exs; 
Eys=vy  1  *I(mode+l :2*mode);  Eysl(m)=Eys; 

etet=-j*k*eta/2/k*(Exs*costet*cosphi+Eys*costet*sinphi);  etetl(m)=etet; 
ephi=-j *k*eta/2/k*(-Exs*sinphi+Eys*cosphi);  ephi  1  (m)=ephi; 

etetrcs(m)=(-j  *k*eta/2/k)*Exs; 
ephircs(m)=(-j  *k*eta/2/k)*(Eys); 

end; 

%%%%%%%%  Calculation  of  RCS  %%%%%%%%% 
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fori=l:mode 

RCStet(i)=10*loglO(abs(etetrcs(i)).A2*2*k+le-10); 
RCSphi(i)=10*loglO(abs(ephircs(i)).A2*2*k+le-10); 
end; 

RCS=reshape(RCStet,modex,modey); 

u=reshape(u,modex,modey);  v=reshape(v,modex,modey); 

mesh(u,v,RCS); 

%  All  the  data  are  saved  for  use  in  synthesis  program. 

save  b24  RCS  u  v  mode  modex  modey  Exsl  Eysl  V  Z  deltax  deltay  xnxmyn  ym  Rs  I 

etetrcs  ephircs  etetl  ephil  etet  ephi  phii  teti  RCSphi 
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%%%%%%%%%%%%%0/o0/o0/o0/o0/o0/o%%0/o%0/o0/o0/o0/o%0/o0/o0/o0/o%%0/o%%%%%% 

%  syntm  % 

%  % 

%  This  program  calculates  the  synthesis  results  for  bistatic  case.  % 

%  The  result  will  be  the  resistivity  which  is  found  % 

%  from  the  synthesis  equations.  % 

%  % 

%  Ugurcan  Samli,  July  1996  % 

%  % 

%%%%%%%%%0/o0/o%0/o%0/o0/o%%0/o%%%0/o0/o%%%0/o0/o0/o0/o%%%%0/o0/o0/o0/o%% 


clear  all; 
format  long; 

%  Reading  in  the  data  file  from  the  scatterring  program  (platedcs.m) 

load  b24.mat; 
disp('Data  file  is  loaded.'); 

eta=377;  lambda=l;  t=pi/180;  k=2*pi; 
uxx=linspace(- 1 , 1  ,mode); 
vyy=linspace(- 1 , 1  ,mode); 

%  Reading  the  E-field  components  from  the  data  file. 

efld(l  :mode)=Exsl(:);  efld(mode+l  :2*mode)=Eysl(:); 
rs=[Rs,Rs];  rs=diag(rs); 

%%%%%%%%%%  Filling  the  matrices  %%%%%%%%%%%%% 

%  Original  resistivity  is  subtracted  from  the  impedance  matix. 

Z=Z-rs; 

form=l:mode 

forn=l:mode 

argy=k*deltay*v(m)/2; 
argx=k*deltax*u(m)/2; 
sincx^l;  sincy=l; 
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if  abs(argx)>le-5,  sincx=sin(argx)/argx;  end; 
if  abs(argy)>le-5,  sincy=sin(argy)/argy;  end; 
expxy=exp(j*k*(xn(n)*u(m)+yn(n)*v(m))); 

fhml  (m,n)=expxy*sincx*sincy*deltax; 
fhm2(m,n)=expxy  *  sincx*  sincy  *  deltay ; 

end; 

end; 

%  Filling  the  scattered  matrix  (First  equation),  and  current 
%  coefficients,  I. 

FFl=pinv(fhml);  FF2=pinv(mm2); 

FFI=[FF  1  ,zeros(mode,mode);zeros(mode,mode),FF2]; 

I=(efld*FFI.'); 

%%%%%%%%%%%  Matrix  Operations  %%%%%%%%%%% 

%  Second  equation  (Equation  63)  operations.  Excitation  vector  V  is  from  the 
%  data  file. 

ZIxy=Z*I.'; 

AA=V.'-ZIxy; 

ZL=(AA./I.'); 

ZLl=real(ZL); 

ZLx=(deltay/deltax).  *ZL  1  ( 1  :mode); 
ZLy=(deltax/deltay).*ZLl  (mode+1 :2*mode); 

%  Reconstructed  resistivity  rx. 

rx=reshape(ZLx,modex,modey ) ; 
rorg=reshape(Rs,modex,modey); 
rx=rx-rorg; 

uxx=reshape(uxx,modey,modex);  vyy=reshape(vyy,modex,modey); 
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%%%%%%%%%%  Plots  %%%%%%%%%%%%%%%%%%%% 

%  Plot  of  residue  (rx)  resistivity. 

mesh(uxx',vyy,rx);  grid;  view(30,30); 

%save  si 2  uxx  vyy  rx  rorg 

title('Comparision  of  the  Resistivities;  in  x-direction'); 

xlabel('u'),  ylabel(V),  zlabel('Resistivity  Differences'); 
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function  [nu,nv]=dcsin(mode) 

%%%%%%%%%%%%%%%%%%0/o%0/o0/o0/o%0/o0/o0/o0/o0/o0/o0/o0/o%0/o%0/o%0/o0/o0/o0/o 

%  This  function  calculates  the  number  of  points  inside  the  % 

%  unit  circle,  and  shows  these  points.  % 

%  % 

%  function  [nu,nv]=dcsin(mode)  % 

%  % 

%  Ugurcan  Samli  % 

%%%%%%%%%%%%%%%0/o0/o%0/o%0/o0/o%0/o0/o%%%0/o%%%0/o%0/o%0/o%%%0/o 

%  Area  ratio  of  the  unit  circle  and  the  total  area. 

tot_pts=l  .28*mode; 

totsqr=sqrt(totjpts); 

if  abs(round(totsqr)-totsqr)<0.2 

nu=round(totsqr); 

nv=round(totsqr); 
else 

nu=ceil(totsqr); 

nv=floor(totsqr); 
end; 

[u,v]=center(nu,nv); 

n_uv=nu*nv; 

sintet=sqrt(u.A2+v.A2); 

fors=l:n_uv 

ifsintet(s)>l 

u(s)=0; 

v(s)=0; 
else 

u(s)=u(s); 

v(s)=v(s); 
end; 

end; 

aa=find(u=0);  bb=find(v=0); 
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u(aa)=D;  v(bb)=n; 

[nu  nv  length(u)  length(v)] 

plot(u,v,,+,),xlabelCu'),ylabel(' v');  axis([- 1,1,-1,1]); 
title(Toints  inside  the  unit  circle'); 
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function  [x,y]=center(modex,modey) 

%%%%%%%%%%%%%%%%%%0/o0/o%%%%%0/o%0/o%0/o0/o0/o%0/o%0/o0/o0/o0/o%0/o 

%  center.m  % 

%  This  Function  finds  the  center  points  of  the  patches  % 

%  in  the  Direction  Cosine  Space.  % 

%  % 
%  [x,y]=center(modex,modey) 
% 

%  % 

%  Ugurcan  Samli,  May,  1 996.  % 

%%%%%%%%%%%%%%%%%%0/o%0/o0/o0/o0/o0/o0/o0/o%0/o%0/o%%%0/o%0/o0/o0/o%0/o 

Lx=2;  Ly=2; 

deltax=Lx/modex;  deltay=Ly/modey; 

a=0; 

for  s=l:modey 
for  ss=l:modex 

a=a+l; 

x(a)=(-Lx/2+deltax/4+(ss- 1  )*deltax); 

y(a)=(-Ly/2+deltay/4+(s-l)*deltay); 

end; 
end; 
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%  green.m 

%  Ugurcan  Samli,  Resistive  Plate  MM  solution 

%  This  function  computes  the  integral  of  green's  function. 

function  gg=green(a,b,c,d,g,h,int); 

w=0.225*h; 

ifint=0 

r=sqrt(((b-a)A2+(c-d)A2)+wA2); 
gg  1  =exp(-2*j  *pi*r)/r; 
gg=ggl*g/4/pi*g; 
end 

ifint=2 

xg=[0.5773503  -0.5773503]; 

ag(l:2)=[l,l]; 
end 

ifint~=0 

int2=int/2; 
for  s=l:int2 

sl=int-s+l; 
ag(sl)=ag(s); 
xg(sl)=xg(s); 
xg(s)=-xg(s); 
end; 

xhigh=a+g/2; 
xlow=a-g/2; 
bb=(xhigh-xlo  w)/2 ; 
aa=(xhigh+xlow)/2 ; 
sum=0; 

for  i=l:int 

x=bb*xg(i)+aa; 
r-sqrt((b-x)A2+(c-d)A2+wA2); 
sum=sum+ag(i)*exp(-j  *2*pi*r)/r; 
end; 
gg=sum*bb/4/pi*g; 
end; 
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%  green  l.m 

%  Ugurcan  Samli,  Resistive  Plate  MM  solution 

%  This  function  computes  the  integral  of  green's  function. 

function  gg=greenl(a,b,c,d,g,n,int); 

w=0.225*g; 

ifint==0 

r=sqrt(((a-b)A2+(c-d)A2)+wA2); 
ggl=exp(-2*j*pi*r)/r; 
gg=ggl*h/4/pi/h; 
end; 

ifint=2 

xg=[0.5773503  -0.5773503]; 

ag(l:2)=[l,l]; 
end; 

ifint~=0 

yhigh=c+h/2; 
ylow=c-h/2; 
bb=(yhigh-y  lo  w)/2 ; 
aa=(yhigh+ylow)/2; 
sum=0; 
fori=l:int 

y=bb*xg(i)+aa; 
r=sqrt((a-b)A2+(y-d)A2+wA2); 
sum=sum+ag(i)*exp(-j  *2*pi*r)/r; 
end; 
gg=sum*  bb/4/pi/h; 
end; 
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%  green2.m 

%  Ugurcan  Samli,  Resistive  Plate  MM  solution 

%  This  function  computes  the  integral  of  green's  function. 

function  gg=green2(a,b,c,d,g,h,int); 

w=0.225*g; 

ifint=0 

r=sqrt(((a-b)A2+(c-d)A2)+wA2); 
ggl=exp(-2*j*pi*r)/r; 
gg=ggl*h/4/pi/h; 
end; 

ifint=2 

xg=[0.5773503  -0.5773503]; 

ag(l:2Hl,l]; 
end; 

ifint~=0 

int2=int/2; 

for  s=l:int2 
sl=int-s+l; 
ag(sl)=ag(s); 
xg(sl)=xg(s); 
xg(s)=-xg(s); 
end; 
yhigh=c+h/2; 
ylow=c-h/2; 
bb=(y  high-y  low)/2 ; 
aa=(yhigh+y  low)/2 ; 
sum=0; 
for  i=l:int 

y=bb*xg(i)+aa; 
r=sqrt((a-b)A2+(y-d)A2+wA2); 
sum=sum+ag(i)*exp(-j  *2*pi*r)/r; 
end; 
gg=sum*bb/4/pi/h; 
end; 
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%  green3.m 

%  Ugurcan  Samli  Resistive  Plate  MM  solution 

%  This  function  computes  the  integral  of  green's  function. 

function  gg=green3(a,b,c,d,g,h,int); 

w=0.225*h; 

ifint=0 

r=sqrt(((a-b)A2+(c-d)A2)+wA2); 

ggl=exp(-2*j*pi*r)/r; 

gg=ggl*g/4/pi/g; 
end; 

if  int=2 

xg=[0.5773503  -0.5773503]; 

ag(l:2)=[l,l]; 
end; 

ifint~=0 

int2=int/2; 
for  s=l:int2 

sl=int-s+l; 
ag(sl)=ag(s); 
xg(sl)=xg(s); 
xg(s)=-xg(s); 
end; 
xhigh=a+g/2; 
xlow=a-g/2; 
bb=(xhigh-xlow)/2; 
aa=(xhigh+xlo  w)/2 ; 
sum=0; 
for  i=l:int 

x=bb*xg(i)+aa; 
r=sqrt((x-b)A2+(c-d)A2+wA2); 
sum=sum+ag(i)*exp(-j*2*pi*r)/r; 
end; 
gg=sum*bb/4/pi/g; 
end; 
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%  green4.m 

%  Ugurcan  Samli  Resistive  Plate  MM  solution 

%  This  function  computes  the  integral  of  green's  function. 

function  gg=green4(a,b,c,d,g,n,int); 

w=0.225*h; 

ifint=0 

r=sqrt(((a-b)A2+(c-d)A2)+wA2); 

ggl=exp(-2*j*pi*r)/r; 

gg=ggl*g/4/pi/g; 
end; 

if  int=2 

xg=[0.5773503  -0.5773503]; 

ag(l:2)=[l,l]; 
end; 

ifint~=0 

int2=int/2; 
for  s=l:int2 

sl=int-s+l; 
ag(sl)=ag(s); 
xg(sl)=xg(s); 
xg(s)=-xg(s); 
end; 
xhigh=a+g/2; 
xlow=a-g/2; 
bb=(xhigh-xlow)/2; 
aa=(xhigh+xlow)/2; 
sum=0; 
for  i=  hint 

x=bb*xg(i)+aa; 
r=sqrt((x-b)A2+(c-d)A2+wA2); 
sum=sum+ag(i)*exp(-j  *2*pi*r)/r; 
end; 
gg=sum*  bb/4/pi/g; 
end; 
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%%%%%%%%%%%%0/o%%%%0/o0/o0/o0/o0/o0/o0/o0/o%0/o0/o%%%0/o0/o%0/o0/o%%0/o%% 

%  % 

%  sample.m                                                                   % 

%  This  program  analyses  the  frequency  content  of  the  E-field                             % 

%  in  u  and  v  directions.                                                                                         % 

%  % 

%  Ugurcan  Samli,  August  96                                                                                % 

%  % 

%%%%%%%%%%%%%%%%%0/o0/o%%0/o0/o0/o0/o0/o%%0/o0/o0/o0/o%%0/o%0/o%%%% 

clear  all; 
load  b24.mat 

Efld=Exsl; 
Efld=reshape(Efld,modex  1  ,modey  1 ) ; 

%  m  :  number  of  FFT  points,  k  :  Spatial  Frequency  index. 

%  us  :  Sampling  interval  in  DCS. 

%  This  is  for  bistatic  case  and  normal  incidence. 

m=128; 

cl=0:fh-l; 

us=2/(modexl-l); 

k=cl*2/(m*us); 

%  fft  in  1 -dimension.  FFT  is  performed  only  in  1  column  (or  row). 

k=k(l:fh/2+l); 

el=Efld(6,:); 

e2=abs(ffl(el,fh)); 

figure(l); 

plot(k,e2(l:fh/2+l));grid; 

xlabelCNumber  of  samples  in  DCS'); 

ylabel('abs(fft)');  title('3  by  3  plate'); 

%  fft  in  2-dimension. 

fNabs(ffi2(Efld,fh,m)); 

[ku,kv]=meshgrid(k,k) ; 

figure(2); 

mesh(ku,kv,ff(l:fh/2+l,l:fh/2+l)),  grid;  view(30,30); 
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xlabel('  Number  of  Samples  in  u-direction'); 
ylabelCNumber  of  Samples  in  v-direction'), 
zlabel('abs(fft)');  title('3  by  3  plate'); 
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